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ABSTRACT 

Ices regulate much of the chemistry during star formation and account for up to 5U7o 
of the available oxygen and carbon. In this paper, we use the Spitzer c2d ice sur- 
vey, complimented with data sets on ices in cloud cores and high-mass protostars, to 
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determine standard ice abundances and to present a coherent picture of the evolu- 
tion of ices during low- and high-mass star formation. The median ice composition 
H20:CO:C02:CH30H:NH3:CH4:XCN is 100:29:29:3:5:5:0.3 and 100:13:13:4:5:2:0.6 to- 
ward low- and high-mass protostars, respectively, and 100:31:38:4:-:-:- in cloud cores. 
In the low-mass sample, the ice abundances with respect to H2O of CH4, NH3, and the 
component of CO2 mixed with H2O typically vary by <25%, indicative of co-formation 
with H2O. In contrast, some CO and CO2 ice components, XCN and CII3OII vary by 
factors 2-10 between the lower and upper quartile. The XCN band correlates with CO, 
consistent with its OCN~ identification. The origin(s) of the different levels of ice abun- 
dance variations are constrained by comparing ice inventories toward different types of 
protostars and background stars, through ice mapping, analysis of cloud-to-cloud vari- 
ations, and ice (anti-)correlations. Based on the analysis, the first ice formation phase 
is driven by hydrogenation of atoms, which results in a H20-dominated ice. At later 
prestellar times, CO freezes out and variations in CO freeze-out levels and the subse- 
quent CO-based chemistry can explain most of the observed ice abundance variations. 
The last important ice evolution stage is thermal and UV processing around protostars, 
resulting in CO desorption, ice segregation and formation of complex organic molecules. 
The distribution of cometary ice abundances are consistent with with the idea that most 
cometary ices have a protostellar origin. 

Subject headings: astrochemistry; stars: formation; ISM: molecules; molecular pro- 
cesses; circumstellar matter; ISM: abundances; ISM: lines and bands; infrared: ISM 



1. Introduction 

Star formation begins with the collapse of an interstellar cloud core to form a protostar, which grows 
from infall of envelope material and later through an accretion disk, where planets may form. This 
physical evolution from cores to planetary systems is accompanied by a chemical evolution, which 
will affect planet and planetesimal compositions. Much of this chemical evolution takes place in 
icy mantles on interstellar grain surfaces, and the aim of this paper is to consolidate the knowledge 
of ice evolution during star formation coming out of the Spitzer Space Telescope mission. 



The first ices were detected in the interstellar medium almost 40 years ago (iGillett fc Forrest 



19731 ) and from previous work, mostly toward high-mass protostars, H2O, CO and CO2 ices are 
known to be common during the cold and dense stages of star formation, with abundances reach- 
ing 10~^ nun- The ice mantles also contain smaller amounts of CII3OH, CH4, NH3 and XCN 



fe.glCibb et al.l 20041) . CH3OH and other ices are proposed sources of complex organic molecules 



(jCharnlev et al.lll992l : iGarrod et al.ll2008l : lOberg et al.ll2009bl ) and determining ice abundances and 



production channels in star forming regions is therefore of great interest for studies of prebiotic 
chemistry. Most ices except for CO are predicted to form in situ on interstellar grain surfaces 



through hydrogenation and oxygenation of atoms and small molecules (jTielens &: HagenI Il982l) 



Heat and UV radiation from protostars may result in additional ice processing and iGibb et al 



(|2004l ) suggested that some CO2 ice, and all CH3OH and XCN ice form through such processing 
based on observations toward high-mass protostars with the Infrared Space Observatory (ISO). This 
has been challenged by observations of abundant CH.cjOH and XCN ice toward low- mass protostars 
(jPontoppidan et al.ll2003bl : Ivan Broekhuizen et alJl2005l ). where ices are protected from strong UV 
fields during most of their life time. 

Spitzer^s high sensitivity made it possible to observe ices t oward low mass protostars and 
background stars at 5-30 urn.. The c2 ( i program (lEvans et al.ll2003l) obtained IRS spectra of 50 low- 



mass protostars (JBoogert et al.ll2008l : iPontoppidan et al.l l2008l : lOberg et al.l 120081 : iBottinelli et al 



2OIOI . from now on Paper I-IV), providing an unprecedented sample size that complements the ISO 
data on high-mass protostars. Additional Spitzer ice observations e xist on ices toward background 



sources, looking through molecular clouds at a range of extinction s (JBergin et al.ll2005l : iKnez et al 



2005 



Pontoppidan 



UV environment (Reach et al. 



2006: 



200! 



J" 



Whittet et al.ll2009l : iBoogert et al. 



2011 



and toward protostars in a high- 



Because the Spitzer spectrometer had a lower cutoff at A ~ 5 /im, complete ice inventories can 
only be obtained by adding complementary ground-based spectroscopy to cover the strong 3 //m 
H2O, the 4.65 /im CO, the 4.6 iim XCN, and the 3.53 fj,m CII3OH features. The latter feature has 
been used to validate the use of the 9.7 /im feature to derive CH3OH ice abundances (Paper I). The 
XCN feature has been shown empirically to consist of two components, w hose relative contributions 
to the feature vary from source to source (Ivan Broekhuizen et al.ll2005l ). One of the components 
can be securely assigned to OCN~ from comparison with laboratory spectra, while the carrier of 
the second component, peaking at 2175 cm^^, is unknown. Suggested carriers of the 2175 cm^^ 
component include OCN~ present in a di fferent ice environment^ resulting in a different spec trum. 



as well as completely different molecules (jPendleton et al. 



1999 



van Broekhuizen et al.ll2004l ). 



In the analysis of astrophysical ice spectra, comparison with laboratory ice spectra is key, both 
to identify ice species (such as OCN~ above) and to characterize the ice morphology; i.e., some ice 
bands, such as CO and CO2, canno t be fitted by a single laboratory ice mixture, but seem t o trace 



molecules in two or more ice phases ([Tielens et al.lll991 



Chiar et al. 



1995l : lDartois et al.lll999l . Paper 



II) . Traditionally the observed spectra have been directly compared to a superpositio n of pure and 



mixed ice spectra, characterizing the ice morp hology in each source independently (e.g. lMerrill et al. 



1976 



Gibb et al.ll2004l : IZasowski et al.ll2009l ). The constraints are often degenerate, however, since 



ice spectral features vary with ice composition, temperature and radiation processing. To address 
this, ice bands ha ve instead been decomposed phenome nologic ally into a small s e t of unique spectral 
compo nents (e.g. iTielens et al.lll99ll : iPendleton et al. 1999; Keane et al.ll200ll : IPontoppidan et al. 
2003al . and §2.2). The contributions of the derived components are then used to characterize 



the spectral bands in all sources in the sample. This method was employed in papers I and II. 
Because all observed spectra are decomposed into the same, small number of components this 
method provides information on the sample as a whole, i.e. it directly shows which parts of the 



spectral profile are ubiquitous and which are environment dependent. This is crucial information 
when assigning a component carrier - without this, the degeneracy is almost always too large to 
draw conclusions about the structure of the ice from spectral profiles alone. 

Building on the analysis in Paper I-IV, this paper aims to establish an ice evolution sce- 
nario that also take s into account the present knowledge of ice abundances toward cloud cores 
([Boogert et al.ll201ll ) and high-mass protostars (jGibb et al.l 12004 ). Section [2] summarizes the sam- 
ple characteristics. Section [3] presents new ice data for 10 low-mass protostars, followed by ice 
abundance medians toward low- and high mass protostars. The variation among the low-mass 
protostellar abundances are investigated followed by comparison of this ice sample with ices toward 
high-mass protostars, low-mass protostars in a high-radiation environment, and background stars. 
The reasons for the abundance variations seen for some ice features are further explored through 
analysis of spatial differences and correlation studies. Further, new data are used to test previous 
hypotheses about the carrier(s) of the XCN band. The results are discussed in ^ with respect 
to different ice formation scenarios and ice chemistry in low-mass versus high-mass star-forming 
regions. 



2. Sample and ice features 
2.1. Sample 



Spitzer-IRS spectra were obtained for ~50 low-mass protostars with ice features as part of the c2d 
Legacy program (PIDs 172 and 179), a dedicated open time program (PID 20604) and a Guaranteed 
Time Observations (GTO) program (PI Houck). All sources were included in Papers II-III on CO2 
and CH4, while the c2d sources alone were investigated in Paper I and IV. 



The Spitzer spectra were complemen ted by ground-based Keck/NIRSPEC (JMcLean et al 



op 
19981 ) and VLT/ISAAC (|Moorwoodl [l997l ) L-band spectra of the H2O ice 3 fim feature. In ad- 
dition, the combined sample partly overlaps with a 3-5 p,m. VLT survey of C O and XCN ice toward 



low-mass protostars (jPontoppidan et al.ll2003bl : lvan Broekhuizen et al.ll2005l ) and additional CO ob- 
servations were obtained with Keck/NIRSPEC (Paper II). The VLT survey contains many sources 
in Ophiuchus, for which the Spitzer data are taken from GTO programs. These sources are key for 
a comprehensive investigation of how the XCN complex, composed of OCN^ and potentially a sec- 
ond carrier, relates to other ices. The reduction of the spectra to obtain abundances are described 
in the individual c2d papers and the additional GTO sources have been reduced using an identical 
procedure to those presented in Papers I and IV. 

The lo w-mass sarnple is complemented with 9 ISO sources, representative of high-mass proto- 
stars, from ICibb et al.l ( 2004 ) which were re-analyzed in Papers I-III to ensure that the low- and 
high-mass ice abundances are consistently derived using the same component analysis and band 
strengths. The sources were not included in Paper IV, and NH3 abundances are therefore taken 



from lGibb et alj (|2004l ). In summary, our sample consists of 63 YSOs, that have been a nalyzed in a 



homo geno us way. When compa red with ice abundances toward background stars from iKnez et al 



(|2005l ) and lBoogert et alJ (|201ll ). this sample spans all evolutionary stages from molecular clouds to 
disks and six orders of magnitude in stellar luminosities and a range of star-forming environments. 



2.2. Ice features 

Figure [T] shows the main ice features seen toward background stars and protostars with the 
identifications of different bands to H2O, CO, CO2, CH3OH, NH3, CH4 and OCN~ ice marked. 
H2O column densities are derived from the 3 fim band whenever such data are available; otherwise 
the librational band (10-30 fxm.) is used. CH3OH and NH3 abundances are generally derived from 
the 9-10 fj,m bands. OCN~ abundances are derived from one of the two components (centered at 
2165 cm~^) that together make up the XCN band (the total band centered at 4.62 /zm next to 
the CO feature). The OCN~ identification is based on laboratory spectra. The abundance of the 
second XCN carrier (associated with the '2175 cm~^' band) is derived assuming the OCN~ band 
strength. 

The CO and CO2 spectral features can be deco mposed into a nurnber of components corre- 
sponding to CO and CO2 in different ice mixtures (jPontoppidan et al.l l2003bl . Paper II). CO is 
decomposed into three components corresponding to pure CO ice, CO mixed with CO2 (C0:C02) 
and CO mixed with II2O ice (C0:Il20). CO2 is similarly decomposed into four components corre- 
sponding to pure CO2 ice, CO2 mixed with CO (C02:C0), CO2 mixed with H^^O ice (CO^rH ^O) 
and a shoulder which has been associated with CO2 mixed with CH3OH (jDartois et al.lll999l ). 



More tentative identifications of spectral bands to HCOOH, NIl|, CII3CH2OH and CII3CIIO 
are also marked in Fig. [TJ Most of these tentative identifications are found in a complex band 
between 5 and 7 /xm. To constrain the carriers of this band it was decomposed into five different 
components (C1-C5) after subtraction of the contribution from II2O ice in Paper I. CI and C2 make 
up the 6 ;um band, C3 and C4 the 6.8 fim band and C5 is a broad, underlying feature that covers 
the entire 5-7 /im region. From comparison with laboratory ice spectra, CI has been identified 
with HCOOH and H2CO, C2 with HCOO" and NH3, C3 with NH+ and CH3OH, C4 with NH+ 
and C5 with warm H2O and anions (Paper I). 



3. Results 



Figure [T] shows the similar 3-20 /xm spectra toward protostars spanning orders of magnitude in 
luminosities, suggestive of a partly shared ice formation history during all types of star formation. 
The spectra are not identical, however, and median abundances of identified ices toward low- and 
high-mass protostars are presented in §3.2. Ice abundance variations toward low-mass protostars are 
further investigated in §3.3 to constrain which ice form together with H2O and which do not. In §3.4 



the low-mass sample is contrasted with ice abundance distributions toward high-mass protostars 
to test the significance of any differences noted in §3.2. §3.5 continues with a comparison of ice 
distributions between clouds and protostars to constrain which ices require protostellar processing 
to form. §3.6. presents an analysis of spatial variations in ice abundances, both on large-scale 
between clouds, and on small scales within a single core. Finally, ice correlations are used in §3.7 
to provide further constraints on XCN formation, and H2O ice anti-correlations are presented to 
investigate which ices form in competition with H2O. 



3.1. New ice abundances in Ophiuchus 

To enlarge the sample of low-mass sources with available XCN data, a set of GTO data in Ophiuchus 
was analyzed in the same way as the c2d sources. The five components, Cl-5, from Paper I are 
fitted to the 5-7 /im complex with their peak optical depths as free parameters. The NH3 and 
CH3OH abundances toward the same sources are determined from their 9.0 and 9.7 /im features, 
using a 4*^ order polynomial to remove the underlying silicate feature (Paper IV). After continuum 
subtraction, column densities are derived by fitting two Gaussians to the observed spectra around 
the expected band positions and using the same band strengths as in Paper IV (Tabled]). The 
flux and optical depth spectra are shown in Appendix A together with fitted peak positions and 
band widths. While NH^ is not explicitly discussed here except in the context of the N-budget, 
maximum abundances for the Ophiuchus sources can be extracted from Table [J using the same 
conversion factor as in Paper I. HCOOH abundances are not reported, since the previously used 
feature may have a significant contribution from other carriers (Appendix B). The presence of a 
feature at 7.25 /Ltm is in itself evidence for the formation of complex ices, even though its dominating 
carrier (probably HCOOH or CH3CH2OH) is speculative. 



3.2. Median abundances 

Figure[2]shows the median ice abundances with respect to H2O ice for low- and high-mass protostars, 
calculated from Paper I-IV and from the new Ophiuchus values reported in Table [H All ice 
abundances are expressed with respect to H2O because H2O ice forms early during star formation, 
is the most abundant ice, and has a high sublimation point. The median values are calculated: 1) 
using only the detections, and 2) using the Kaplan-Meier (KM) estimate of the survival function. 
The latter is a non-parametric procedure that takes into account the constraints provided by upper 
limits. When calculating the KM estimate the detections and upper limits are ordered from low 
to high and the upper limits are given the values of the nearest lower detections. For example 
in a sample of four detections of 0.5, 1, 3 and 4 and an upper limit of <2, the upper limit is 
treated as a detection of 1 . The KM estimate and how to apply statistical tests on it is reviewed by 



Feigelson &: NelsonI (119851 ). The medians calculated from the KM estimate provide more accurate 



descriptions of ice populations with significant upper limits as demonstrated in Fig. [21 where the 



low- and high-mass CH3OH medians and the NH3 and CH4 median toward high- mass protostars 
are 40-70% lower when taking into account the upper limits. 

Table [2] lists the medians that best describe the data set for total ice abundances, while Table 
[3] provides medians for all ice components together with their lower and upper quartile values. 
Quartiles are favored over standard deviations since most ice abundances do not follow a normal 
distribution (see below). For the species where the Kaplan- Meier estimate is >20% lower than 
the detection-based median, both are listed. CO and CO2 are the most abundant ices after H2O, 
with low-mass protostellar median abundances of 29% with respect to H2O. NH3, CH3OH and 
CH4 have comparable median abundances of 3-5%, while XCN is rare at <1%. The CH3OH and 
NH3 abundance medians are similar toward low- and high-mass protostars, when the upper limits 
toward high-mass protostars are taken into account, while CO, CO2 and CH4 are more abundant 
and XCN is less abundant toward low-mass protostars. The low CO2 abundances in this high-mass 
sample are consistent with recent observations of CO2 ice abundances of 10-18% to ward high-mass 



protostars in the g alactic center, but not toward high-mass YSOs in the LMC ()An et al. 



200s 



Oliveira et al.ll2011l ). Testing whether perceived differences between low- and high-mass protostars 
are significant requires an investigation of the ice abundance distributions within each sample of 
protostars, however, and this is done in §3.4 after presenting ice abundance variations toward 
low-mass protostars in §3.3. 



3.3. Ice abundance variations 

Ice abundance variations between different sources depend on how sensitive the ice formation and 
destruction pathways are to the local environment. Investigating protostellar abundance variations 
thus puts constraints on when and where different ices form (variations are not dominated by 
reported abundance uncertainties from fitting the observed spectra). Figure [3] demonstrates that 
some ice abundance distributions are skewed with a high-abundance tail. All abundance histograms 
are therefore log-transformed and then centered on the median detected low-mass protostellar ice 
abundance, with bin sizes proportional to the low-mass abundance variances. For ices where the 
median is unaffected or barely affected (<20%) by including upper limits, the detections alone are 
shown. For ice abundances that are better constrained by including upper limits (CO, CH3OH, 
XCN, OCN^, 2175 cm-\ C2 and C5 toward low-mass protostars and CH3OH, NH3, CH4, XCN, 
C0:C02 and C5 toward high-mass protostars), the distributions include 3-0" upper limits. Ap- 
pendix C contains figures that compare the distributions with and without upper limits included. 
Differences between low- and high-mass protostellar abundances are discussed further below, while 
this section focuses on the larger low-mass sample. 

Since ice abundances are with respect to H2O, a small variation of a species indicates co- 
formation with H2O, while large ice abundance variations are indicative of different formation 
and/or destruction dependencies than H2O ice. Figure H] shows that the total CO2, CH4 and 
NH3 abundance distributions are narrow, while CO, OCN~, the 2175 cm~^ XCN component and 



CH3OH have broader distributions with order of magnitude abundance variations between different 
sources. Figure [5] shows that the narrow CO2 distribution is due to C02:H20 ice, while all other 
CO2 and CO component distributions are broad. The pure CO and CO2 ice components have 
the broadest distributions, consistent with their expected dependence on the protostellar envelope 
temperature for ice evaporation and segregation. The 5-7 /im components, Cl-5, span the full 
range of variations observed among the identified ice components, including the previously noted 
narrow C3 distribution (Paper I). 

Protostellar ice heating was explored in Paper I and II as explanations for observed abun- 
dance variations, using pure C0/C0:H20 and pure C02/C02:H20 as diagnostics of the amount 
of moderate ice heating, and H20/silicate as a probe of complete ice sublimation. Correlations 
be tween H2O /silicate a nd the different Cl-5 components are extensively discussed in Paper I and 



by lBoogert et al.l (|201ll ) for protostars and background stars. The main conclusion was that C3/C4 
does depend on ice sublimation, supporting its identification with low- and high-temperature NH^ 
salts. 

Moderate ice heating is predicted to reduce the abundances of volatile ices (decreasing the pure 
CO abundance), segregate previous ice mixtures (resulting in pure CO2), drive ice crystallization, 
activate acid-base chemistry and cause diffusion of ice radicals and thus the formation of more 
complex species. None of the ice abundances in this study are however correlated with either ice 
temperature tracer (not shown) except for the previously found C5/H2O correlation with pure 
C0/C0:H20 (Paper I). Moderate (transient) ice heating alone seems to play a minor role in simple 
ice formation. This result does not exclude that the combined differences in temperature and UV 
radiation around low- and high-mass protostars results in different ice chemistries. Whether this is 
the case is the topic of next section. 



3.4. Low- versus high-mass protostars 

Figures HHS] show that many ices are similarly distributed toward low- and high-mass protostars, 
but as noted in §3.2, the total CO2, CO and CH4 abundances are lower toward the high-mass 
sources. Figure [5] reveals that the difference in CO2 abundances between low-mass and high-mass 
stars is due to a difference in the C02:H20 ice component, which typically dominate the spectra 
toward low-mass sources (Paper II). All CO components are centered around a lower value toward 
high-mass protostars compared to low-mass protostars. Among the Cl-5 components, C4 seems 
more abundant toward high-mass protostars. 

The significance of the visual differences in Figs. I2|4f [6l and other proposed differences be- 
tween low-mass and high-mass protostars can be evaluated statistically. Kruskal-Wallis one-way 
analysis of variance by ranks is a non-parametric met hod to test the equality of medians in two or 



more groups of observations (JKruskal fc Wallislll952l ). The test statistic does not assume normal 
distributions and is given by 



Sf=l 



riiin - r 



i2 



i^ = (iV-l) /=V\1 ^2' (1) 

where g is the number of groups, rii is the number of observations (ice abundances) in group i (low 
or high-mass protostars), Tjj is the rank of observation j in group i (here the rank of a source in 
terms of the investigated ice abundance), N is the total number of observations across all groups, 
fj is the mean rank of group i and f is the mean rank in the whole population. The null hypothesis 
of equal medians is rejected when K > Xa-n~i where a is the significance level and g — 1 the degrees 
of freedom for the x^ statistic. The Xaq~i values can be looked up in tables and are also part of 
most statistic packages such as R. 

Applying the test to the low-mass versus high-mass detected ice abundances shows that there 
is a significant difference (at the 95% level) between the two samples for CO2, CO, C02:H20, pure 
CO and C4 ice abundances. Differences in C5, CH3OH or OCN~, which have been suggested to 
be more efficiently formed toward high-mass protostars, are not statistically significant. To check 
whether including upper limits affect these results, a similar test (the logrank test) was run in R 
using the previously calculated Keplar-Meier estimates and the function survdijj. The results are 
consistent with the Kruskal-Wallis test, except that CH4 is different in low -and high-mass stars at 
the 98% confidence level when including upper limits. 

In the above analysis, all low-mass protostars are situated in star-forming regions without any 
nearby massive stars. It is therefore difficult to disentangle which differences are intrinsic to the 
low- versus high-mass star-forming processes and which differences depend on the local radiation 
environment. IC 1396A, known as the Elephant Trunk Nebula, is a dense globule, excited by 
the 4 Myr old 06 star HD 20626 7. Ice observation s - H2O, CO2 and Cl-5 - toward four low- 



mass protostars are presented bv iReach et al.l (|2009l ). The C3/C4 ratio toward these protostars 



resembled the ratio found toward high-mass protostars, suggesting that the ice evolution around 
low-mass protostars depends on the external radiation environment. Figure [7] shows a histogram 
comparison for the CO2 and Cl-5 abundances toward the low-mass, high-mass and IC 1396A 
samples. Visually the IC 1396A C3 and C4 abundances seem to better overlap with the high-mass 
sample, while the IC 1396A C5 upper limits and CO2 detections are more consistent with the 
low-mass sample. None of these differences are, however, significant, nor is the difference in C3/C4 
ratio, likely due to the small IC 1396A sample size. 



3.5. Protostars versus background stars 

Comparisons between protostars and background stars provide direct limits on which species form 
before the protostar turns on and starts to heat and irradiate its surroundings. Tables [2] and [3] 



10 



lists the medianqj of species detected toward background stars in Serpens and in individual cores 
(jKnez et al.ll2005l : iBoogert et al.ll201lh . For the isolated core sample we only use those targets with 
well known H2O column densities from 3 ^m observations. The upper limits on CH4, NH3 and 
OCN~ are generally too high to be useful and there are only a handful of CO ice observations. 
Cl-5, CO2 and CH3OH are, however, detected in a larger number of sources and can be used to 
compare cloud and protostellar ices. 

Previous ice studies toward T aurus revealed a low er CO2 ice fraction in clouds compared to 
in low- mass protostellar envelopes (jWhittet et al.l 120071 . Paper II) which together with laboratory 
studies have enc ouraged hypotheses of CO2 formation from energetic processing of CO ice around 



protostars (e.g., Iloppolo et al.ll2009l). Similarlv. CH3OH ice has been propo sed to be a product o f 



protostellar UV ice-processing (JGibb et al.ll2004l ). Figure El using data from lBoogert et al.l (J2011 



shows that there is no evidence for different CO2 and CH3OII ice abundances toward protostars and 
background stars. Applying the Kruski-Wallis and Keplar-Meier based tests confirm this lack of a 
difference in ice composition between our protostellar ar id cloud sources. Fig ure [8] also shows the 
clear difference between Taurus CO2 ice abundances from lWhittet et al.l (120071 ) and our background 
star ice sample. From this difference, ice abundances toward Taurus background stars are not good 
templates for ice abundances toward average cloud cores. 

In summary, observable ices, except for pure CO2 ice and the C5 component, are present at 
the same abundances in the pre-stellar and protostellar phase. Most protostellar ice abundance 
variations reported in §3.3 must then be due to differences between clouds and to different pre-stellar 
ice formation processes, rather than thermal and UV processing by the protostar. 



3.6. Spatial differences within the lovi^-niass protostellar sample 

3.6.1. Cloud-to- cloud variations 



Variations in ice formation and destruction efficiencies in the pre-stellar stage may be associated 
with different cloud structures and initial chemical conditions in different star-forming regions. This 
is investigated by dividing the low-mass protostellar sample into six groups with sources belonging 
to the star forming associations in Ophiuchus, Serpens, Corona Australis, Perseus, Taurus and a 
collection of smaller associations. Kruskal-Wallis tests show that when considered as a whole, the 
ice abundance medians are statistically i ndistinguisha ble between these groups. This is consistent 
with a Principal Component Analysis in lObergj (120091 ) . which revealed no differences between the 
different clouds when simultaneously considering all ice abundances. 

Cloud-to-cloud variations can also be tested for specific ice components. CO2 ice abundances 



^The background star medians are based on detections alone, since applying survival analysis shows that including 
upper limits does not constrain the medians further. 
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in clouds were shown above to be unusually low toward Taurus, while too few CO measurements 
toward background sources outside of Taurus exist to determine whether CO is lower than expected 
as well. Testing the CO and CO2 protostellar differences between the clouds may reveal whether 
prestellar differences are carried over ii ito the protostellar phase. CH3OH is another ice that has 



been observed to be spatially variable (jPontoppidan et al.ll2004l : iBoogert et al.ll201ll ). Differences 



in the pure CO ice component and the CO2 shoulder (often assigned to C02:CH30H interactions) 
between the clouds are indeed statistically significant at a >98% confidence level, while all other 
CO and CO2 components and CH3OH are not. The significant differences exist between the 
low abundances toward Taurus and the smaller associations on one hand, and the more 'normal' 
abundances toward the other clouds (Table H]) . 

The difference between Taurus and the other clou ds is important since much of the current ice 
paradigm is based on analysis of Taurus sources (e.g. IWhittet et al.l 120071 ). The low pure CO ice 
abundances may also explain the low CH3OH ice upper limits toward Ta urus, since CO freez e-out 
is a pre-requisite for CH3OH ice formation without energetic radiation (jCuppen et al.ll2009l ) and 
a low pure CO ice abundance may indicate that catastrophic CO freeze-out has not taken place. 
In summary, some CO ice abundance differences seem to be explained by large-scale differences in 
cloud environments, but overall most cloud and pre-stellar ice formation variability must be due to 
more local effects than cloud-to-cloud variation. Both from cloud and protostellar ice abundances, 
Taurus stands out and care should hence be taken in using Taurus trends as a basis for general 
conclusions on ice evolution during low-mass star formation. 



3.6.2. An ice map toward the Oph-F core 

Within single cloud cores, ice maps may reveal production and destruction pathways that depend 
on the local environment. An ice map of protostars in the Oph-F core has been used previously to 
show that the total CO a nd the C0:H20 abundances decrease monotonically away from the central 
core (jPontoppidanI |2006| ) . The lines of sight probe primarily the dense quiescent core, rather than 
ice in the protostellar envelopes, and the spatial trends were interpreted as catastrophic freeze-out 
of CO in the pre-stellar stage once a certain density and temperature is reached. Figure [9] shows 
that the order of magnitude increase in CO ice with respect to H2O toward the core is accompanied 
by an increase in C02:C0. In contrast, the C02:Il20 ice component is almost constant across the 
core. The XCN band (the 2175 cm~^ feature - OCN~ is not detected) is the only other species that 
increases monotonically toward the densest part of Oph-F. C02:C0 and the 2175 cm~^ band thus 
appear directly related to CO freeze-out. No trend with distance to the core center is seen for any 
other ice component and is also not expected for ices forming early during cloud formation (e.g. 
CII4 and NH3), which are independent of cloud core time scales and CO freeze-out, or of species 
dependent on protostellar heating (e.g. pure CO2 ice) or of components with potentially multiple 
carriers, such as the Cl-5 bands. CH3OH is only detected toward one of the sources and no trend 
can thus be extracted. Still the Oph-F map suggests that CO freeze-out followed by a CO driven 



12 



ice chemistry may account for large pre-stellar ice variations. 



3.7. Protostellar correlation studies 

The trends found for the Oph-F core are here explored for the entire protostellar sample through 
correlation studies between especially CO-chemistry products and H2O ice column densities, and 
between XCN, and CO and CO2 ice components. 



3.1.1. Dependencies on H2O ice column 

In individual cores, the H2O column density has bee n shown to correlate well with dust extinction , 
assumed to trace the total dust mass and H2 column (jWhittet et al.lll988l : lPontoppidan et al.ll2004l ). 
This implies that the average abundance of H2O ice per dust grain is constant across the core. In an 
arbitrary line of sight, the total H2O ice column density depends on both the column density of dust 
grains and on the average number of monolayers of II2O ice on each dust grain. If variations in the 
H2O column density across the sample are dominated by the line of sight fraction of high-density 
material, ices dependent on CO freeze-out are expected to correlate with II2O column density, since 
CO freezes out catastrophically at high densities (§3.6.2). If instead variations in the II2O column 
density across the sample is dominated by the H2O abundance on each grain, the main observable 
trend should be that ices that form in competition with II2O anti-correlate with the II2O column 
density. 

In theory, the dust column density can be estimated from the optical depth of the 9.7 fiui 
silicate feature or from the continuum extinction. Because of a number of practical complications 
we do not attempt to normalize the H2O column density to the dust column, however. First, 
toward protostars, the silicate absorption can be filled in by an unknown amount of silicate emission 
from the warmer parts of the envelope. Second, continuum extinction estimates are affected by 
uncertainties in the intrinsic protostellar SEDs. Third, while the 9.7 /im band correlates well 
with the near-infrared color excess, t he relations vary in different environments by factors of 2 
(|Chiar et al.ll2007l : iBoogert et al.ll201ll ). poss ibly because of variations in the grain size distribution 
or even variations in the grain composition (jvan Breemen et al.l 120111 ). And fifth, the contribution 
of unrelated, ice-less foreground dust is often unknown. 

There are no significant positive correlations (at the 95% confidence level) between HoO and 
any ice abundance, X/H2O, when applying Spearman's rank correlation test (ISpearmanI Il904l ) 
to the low-mass protostellar sample; rank correlation tests are more robust to outliers compared 
to most parametric correlation tests and make no assumptions about whether the correlation is 
linear or non-linear. Figure [TOJ shows that instead, six ice components - CII4, C02:C0, the CO2 
shoulder, C0:Il20, C0:C02 and C4 - are inversely correlated with the II2O column density at 
the 97-99% confidence level. C4 has been associated with NHi" salts and its anti-correlation with 
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H2O may be explained by an increasing salt fraction with ice sublimation because salts have higher 
sublimation temperatures than H2O ice (Paper I). This scenario cannot, however, explain the other 
anti-correlations, since CH4, CO and CO2 are all more volatile than H2O. 

The CO-related anti-correlations may be a result of competitive formation of H2O and CO; 
the more oxygen that is bound up in H2O ice, the less may be available to form gaseous CO 
and thus ices that depend on CO freeze-out. Still in each individual core CO ice abundances and 
H2O ice column densities are expected to correlate because of increasing CO freeze-out toward the 
densest part of the cloud core. The latter may explain some of the scatter in the plots. Finally, 
the CH4-H2O trend is mainly visible for very low H2O column densities (many in Ophiuchus) and 
this relation may be due to more efficient CH4 formation during only the earliest H2O formation 
stage when a large fraction of carbon is still in atomic form. 



3.1.2. XCN correlations 

As stated in §1, the XCN band is composed of two components, peaking at 2165 cm^^ (OCN^) and 
at 2175 cm~^. Toward the Oph-F core, the 2175 cm~^ component is spatially correlated with CO, 
C02:C0 and C0:H20. These correlations were explored in the whole low-mass sample and Fig. 
[TT] shows that the 2175 cm^^ component is significantly correlated with CO and C02:C0 (95% 
confidence with Spearman's rank correlation testjj and the correlation is even stronger (>99%) 
when the high-mass sources are added. The complete XCN band is not correlated with either of 
these components, but it is significantly correlated with the less volatile C0:H20 component (95% 



level). Panel e, finally confirms previous claims in Ivan Broekhuizen et al.l (120041 ) that the relative 



importance of the OCN feature and the 2175 cm ^ feature depend on ice processing, measured 
here by the pure CO / C0:H20 ice ratio. 

These correlations are indicative of a CO-related, single carrier of the entire XCN band with 
a spectral profile that depends on the environment. The simplest explanation is that OCN~ is 
responsible for the entire XCN band (i.e. both the 'OCN~ component' and the 2175 cm^^ com- 
ponent) and that the two components are due to OCN~ in a volatile (CO-rich) ice toward some 
sources, and in a different ice m ixture in warm environments. OCN~ can form at 15 K from HNCO 



in the presence of strong bases (JRaunier et al.ll2003l : lvan Broekhuizen et al.ll2004l ). Consistent with 
the 2175 cm~^ identification with OCN~, this component seems to be correlated (95% level) with 
C4 (one of the two bands ascribed to the base NH^ in Paper I) when only considering low-mass 
protostellar detections. The correlation is not significant, however, when including the high-mass 
sample or the significant upper limits and must therefore be considered as tentative. 



^The Spearman's rank correlation is not yet implemented for censored data in R and the test was therefore not 
repeated including the constraints of the upper limit. Figure [TT] shows such an inclusion would only strongly affect 
the C4 correlation. 
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Discussion 



4.1. C-, 0-, and N-budget 

The amount of C, O and N that are typically bound up in ice mantles is important for the 
life cycle of the elements in the interstellar medium. The total (refractory + volatile) C, O and N 
abundances in the sola r neighborhood are 2.1, 5.8 and 0.58x10"^ per hydrogen nucleus, respectively 
([Przvbilla et al.l 120081 ). The fractional C, O and N abundance in ices with respect to these total 
C, N and O abundances are calculated from the presented ice abundances with respect to H^O ice 



togeth er w ith literature values of the H2O abundance with respect to hydrogen. iPontoppidan et al. 



(|2004l ) and lBoogert &: Ehrenfreundl (J2004l ) modeled the H2O abundance toward four of the low-mass 
protostars in the sample and found riHaO = 4.9ib0.7x lO^^nfj. Toward high-mass protostars the H2O 
abundance can be estimated using the rela tions Ay/rgy = 9 and Nn/Ay ~ 1.9 x 10^^ cm~^ mag~^ 



(Roche Sz Aitken 



1985 



Bohlin et al.lll978l ). It is important to note that this method only provides 



a crude estimate since the silicate feature is affected both by emission and gra in growth; a, recen t 
grain model result in a conver s ion fa ctor that is 30% lower for dense clouds (JEvans et al.ll2009l ). 
Using the values in lcibb et all j^OoJ) for W33A, NGC7538 IRS 9, Mon R2 IRS 3 and S140 IRS 1 



(covering the range of observed ice processing) results in a II2O abundance of (5.0 it 1.9) x 10 ^rin. 



A mean II2O abundance of 5 x 10 
high-mass protostars. 



^riii (i.e. 10 ^ riHa) is therefore assumed for both low and 



The O abundance in the ice is calculated from (xH20+2^co + 2xxco2+^CH30H+a;xcN) X''^H205 
the C abundance from (xco + 2:002 + ^CHsOH + xqu^ + a^xcN) x 'T'H20 and the N abundance from 
(^^NHa + 2;xcN + ^NH+) ^ ''^H2 0) where Xy is the abundance of y with respect to II2O. The median 
upper limit on NH^ is calculated from the total optical depth of the C3 and C4 components, using 
the FW HM of 0.195 and 0.292 /iin , respectively from Paper I and a band strength of 4.4 x 10~^^ 
cm from lSchutte &: Khannal (J2003l ). 



The median percentages of C, O and N bound up in ices, calculated from the ice abundances 
in Table 2, are listed in Table [5j On average 15% of the total (refractory -|- volatile) C, 16% of the 
O and 10% of the N reside in ices toward low-mass protostars. Toward high-mass protostars the 
same calculation yields 12% of the O, 8% of the C and up to 12% of the N atoms. The differences in 
the amount of C and O are due to the low CO and CO2 abundances toward high-mass protostars. 

Table [5] also lists the amount of volatile or non-refractory amounts of C and O that are in 
ices. In the solar neighborhood the non-refractory O abundance is 3.2 x 10~^ per Nh (JMever et al 



19981 ) because a large fractio n of O is bound up in silicates and is therefore not available for ice 
formation (e.g IWhittetll20ld ). Out of this available O, 34% is bound up in ices toward low-mass 
protostars, and 25% toward high-mass stars. The fraction of C in refracto ry material is less well 
const rained, but assuming a non-refractory C abundance of ~ 1 x 10~^ ( Weingartner Sz Draine 
2OOII ) results in 27% of the available C found in ices toward low-mass protostars and 14% toward 



high-mass protostars. 
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The above numbers all relate to the median ice abundances in the sample. The extreme 
ice sources may, however, provide more information on whether ice formation is limited by the 
elemental abundances. IRS 51 and NGC 7538 IRS 9 are calculated to be the most ice-rich sources 
in the sample, when still assuming a II2O ice abundance of 5 x 10~^ with respect to tt-h- Especially 
in the case of IRS 51, the amount of O and C found in the ices approaches the abundances of 
non-refractory O and C (Table [5]). 

It is important to note that additional O and N may be hidden in the ice in the form of 
O2 and N2. Direct constraints on O2 ice abundances are high, 50-100% with respect to CO ice , 
and constraints from O2 interactions with CO ice ar e of a similar order (JEhrenfreund et al.l 11992 : 
Vandenbussche et al.lll999l : lPontoppidan et al.ll2003bl ). O2 ice is, however, unlikely to be abundant 
from the low upper limits on O2 in the gas-phase; it is orders of magnitude less abundant than CO 
and since the two molecules have similar freeze-out and desorp tion properties it is uii l ikely that O2 



is differently partitioned between ice and gas compared to CO ([Goldsmith et al.ll2000l : iPagani et al 



2003 



Fuchs et al.ll200fil ). 



4.2. Ice evolution: early and late pre-stellar ice formation 

Toward low- and high-mass protostars alike, the ices are dominated by H2O, followed by CO2 
and CO (§3.2). C and O are expected to react similarly with hydrogen. The low abundance of CH4 
therefore implies that most C is not in atomic form when the bulk of the H2O ice is forming from 
O early in cloud formation or at low extinctions (Region 1. in Fig. [T2|) . A fast conversion of C into 
CO is also consistent with the CII4/H2O versus H2O anti-correlation toward lines of sight with low 
H2O column densities (§3.7); it suggests that CII4 forms more efficiently in the very beginning of 
H2O formation compared to when the bulk of H2O ice forms at slightly higher extinctions. The 
bulk of H 2O formation is instead associated with continuous formation of CO2, probably from 
CO+OH (jOba et al.ll2010l : lloppolo et al.ll201ll ): the C02:H20 abundance barely varies from source 
to source in the low-mass protostellar sample (§3.3). II2O and CO2 ice mapping confirms this 
scenario ([Bergin et al.ll2005l ). as do the similar CO2 abundances toward cloud cores and protostars 
(§3.5). Most CO that is frozen out in this phase must be converted into CO2. 

From its small abundance variations, NII3 is also inferred to form early with II2O (§3.3). 
The C3 component may be associated with this stage as well and its identification with NH^ 
suggests that there exists a large amount of strong acids that can protonate NII3 in cold ices 



([van Broekhuizen et al.l |2004| ). Anions are notoriously difficult to observe in the ice, however. 



and thus gas phase observations of non-thermally evaporated ices at the edges of clouds probably 
offer the best opportunity to confi rm their existence . It is promising that HCOOH has been 
detected toward translucent clouds ([Turner et al.lll999l ) and a more comprehensive study focusing 
on gas-phase HCOOH and other acids toward cloud edges, where ice formation has began and 
photodesorption is efficient, is warranted. 
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At some point the CO/0 ratio becomes large enough that most frozen-out CO is no longer 
converted into CO2 and H2O is no longer the most efficiently formed ice. We mark this as the 
breaking point between early H2O dominated and late CO-driven pre-stellar ice formation. The 
ices that form in this later stage (Region 2 in Fig. [T2]) are not expected to co-vary with H2O, 
but rather depend on the CO-freeze-out efficiency. The latter is both temperature and density 
dependent and may vary greatly dependent on cloud structure, cloud collapse time scales and 
local radiation environment. This explains why many identified ice and ice components that are 
definitely present during the pre-stellar stage (e.g. CH3OH) have abundances that vary by orders 
of magnitude with respect to H2O (§3.3,3.5) both locally in a single core and when comparing 
different cloud complexes. In addition, the anti-correlation between H2O ice and the abundance of 
several CO-ice components suggest that the time available for this CO-based chemistry decreases 
when a longer time is spent in the H20-ice formation stage, where O is converted into H2O and 
C02:H20 ice (§3.7). 

C02:C0, C0:H20 and the OCN~ feature are all associated with this CO-ice chemistry (§3.5). 
C02:C0 correlates with CO freeze-out in the Oph-F ice map, demonstrating a second, later produc- 
tion pathway of CO2. CO2 may still form from OH-I-CO, but at later stages CO is more abundant 
than O, resulting in s parse H2O formation compared to CO and thus in a CO-rich CO2 ice (see 



Garrod &: Paulv 



2011 



for a similar discussion from a theoretical point of view) . Prom the correla- 
tion studies in §3.7, OCN~ seems to similarly require CO freeze-out and may form from CO-I-NH, 
followed by proton transfer. 



CH3OH and H2CO can form from hydrogenation of CO (ICuppen et al.ll2009l ) and should also 
depend on CO freeze-out levels. There is no strong correlation between CH3OH and CO, however. 
The CO to CH3OH conversion efficiency thus vary in different lines of s ight, probably depe ndent 
on the H/CO r atio, which niay va ry with density and collapse time (jCuppen et al.l |2009| ) and 



cosmic ray flux (IBoogert et al. 



2011 



Ices are continuously exposed to cosmic rays and cosmic-ray 



induced UV photons (jShen et al.l 12004 ) and other 'late' ice features may also depend on cosmic- 
rays, directly or indirectly. In general, observed ice features can however be explained by an early 
ice chemistry phase characteristic of hydrogenation of atoms, followed by atom- addition reactions 
in a CO dominated ice phase. 



4.3. Ice evolution: the protostellar stage 



Though much of the observed ice evolution can be explained by prestellar processes, some ice 
features are formed or destroyed by the protostar. The broadest ice distributions belong to pure 
CO2 and pure CO ice (§3.3), whose variations are best explained by protostellar heating evaporating 
CO ice and distilling or segregating C02-containing ices (Paper Il lOberg et al.ll2009al ). This last 
and most source dependent ice chemistry phase is also pr edicted to result in c omplex ice formatio n 



that requires ice heating for diffusion of large radicals (iGarrod et al. 



2008 



Oberget al.ll2009bl ). 



While the spectra observed to date do not allow for definitive assignments to any specific complex 
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species, the presence of e.g. a feature at 7.25 fj.ni belonging either to HCOOH or CH3CH2OH, 
demonstrates that some complex organic ices are forming efficiently in the pre- or protostellar 
stage. 

Most of the observed differences between low and high-mass protostars can be explained by 
protostellar ice heating (§3.4). CO is comparatively volatile and its low abundances toward high- 
mass sources is consistent with their higher envelope temperatures. The C0:H20 contents toward 
low and high mass objects are comparable, consistent with experiments that show that 5-10% of 
CO can be trapped in H9O ice up to the H2O evaporation temperature (jSandford &: Allamandola 
19901 : ICoUings et aPbooi IPavolle et allboilh . 



The low C02:H20 abundances toward high-mass protostars can also be explained by ther- 
mal ice processing. This component starts to segregate above 30 K resulting in a pure CO2 ice 
component, which desorbs above 50 K. The luminosity of the star will determine the radii for ice 
segregation and desorption and the extent of the 30-50 K shell where pure CO2 ice will be common. 
Toward low-luminosity stars the segregation and desorption radii will be closer to the star resulting 
in more unprocessed ice along the line of sight and thus more C02:H20 as compared to high-mass 
YSOs. 

An alternative explanation for the low C02:H20 fraction toward high-mass stars is different 
cloud formation time scales and t empe rature structures during low and high mass star formation. 
Recent work bv lGarrod &: Paulvl (J201ll ) shows that the switchover between C02:H20 and C02:C0 
ice formation depends on temperature and different cloud temperature structures toward low- and 
high-mass forming clouds may account for the observed differences. The high CO2 abundances 
toward the small sample of low-mass protostars in a high-radiation region (IC 1396A) seems more 
consistent with the protostellar processing scenario, however. 

In the low-mass protostellar sample there is no correlation between pure CO2 ice abundance 
and the total envelope mass as traced by H2O ice (§3.7.1). There is thus no evidence that variations 
in pure CO2 fractions can be explained by variab le envelope sizes (towar d high-mass protostars, 
envelope temperatures and sizes anti-correlate, see Ivan der Tak et al.ll2000l ). Rather large amounts 
of pure CO2 toward some low-mass protostars may be a remnant of past luminosity outbursts 
that heated up the envelope beyond its low-accretion temperature, since ice segregation is an 
irreversible process ( Kim et al.ll2011al )bl). Episodic accretion resulting in luminosity outbursts has 
been invoked to explain the luminosity problem of low-mass YSOs, i.e. the observation that most 
YSOs are significantly less luminous than predicte d by protostellar evolutionary models with steady 



accretion (JEvans et al. 



2009 



Dunham et al.ll2010l ). 



There is no evidence for a higher CH3OH or OCN~ content toward high-mass stars, except 
for the cases of over abundant OCN^ toward W33A and very abundant CH3OH toward GL7009S 
(JDartois et al.lll999l ): thus neither species seems to require stellar UV irradiation or thermal pro- 
cessing to form at typical abundances. The composition of the XCN band varies on average between 
low- and high-mass objects, but as discussed above the carrier is probably dominated by OCN 



m 
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both cases. Of the remaining ice components C4 and C5 seem to have higher abundances toward 
the high- mass protostars, indicative of at least a partial formation pathway involving heat or UV 
irradiation or both, in agreement with the analysis in Paper I. Still, most ice abundances are re- 
markably similar around low- and high-mass protostars, high-lighting the importance of the cold, 
protected stages for ice formation up to the complexity of CH3OH. This is also the conclusion in 
Boogert et all (120 111 ). 



4.4. Ice evolution in disks? 



The protostellar stage includes the formation of an accretion disk, which is the formation site 
of planetesimals and eventually planets. These disks are observed to have cold midplanes, where 
ices should be abundant. It is often difficult to st udy these disk ices direct ly because of confusion 
between disk ices and ices in foreground clouds ([Pontoppidan et al.ll2005l ). Comets in our own 
solar system should, however, carry a record of the ice composition in one such disk midplane. It is 
unknown whether cometary ices reformed in the disk or whether they are at least partly pre-solar or 
protosolar; recent mod els suggest that mu ch of the original ice will in fact survive accretion from the 
envelope into the disk (jVisser et al.ll2009l ). Even if cometary ices originate from protosolar ices, the 
protostellar and protosolar ice abundance populations are not expected to overlap perfectly, since 
the comets originate around a single star, whereas the spread in protostellar ices show the differences 
in ice abundances between different protostars. Assuming that the protosun resembled one of the 
observed protostars, what can be addressed is the question of whether comet ice abundances are 
consistent with a protosolar origin or whether additional ice chemistry in the protosolar nebula is 
required to account for the observed abundances. In light of this. Figure [13] shows a comparison 
betwe en cometary and low-mass protostellar CH .cjOH, CH4 and CO abundances with respect to 
H2O (JBiver et al.l 120021 : iDisanti &: Mummall2008l . and references therein). 



For all three ices, the cometary abundances are well below the median protostellar abundances. 
If the observed cometary ices are pre-solar, this comparison implies either that the ice composition 
in the protosolar envelope was similar to the most carbon-poor ices observed in the protostellar 
sample, or that some ice abundances with respect to II2O were reduced in the solar nebula. Possible 
selective ice destruction processes include desorption of the most volatile ices and UV induced 
chemistry. 

The spread in cometary ice abundances is comparable to the spread found between different 
protostellar sources, suggestive either of a primordial variation in ice abundances within the proto- 
solar envelope or of ice processing in the solar nebula, or a combination of the two. Once a sample of 
cometary CO2 ice abundances become available, the dominant source of variation should be possible 
to constrain, since CO2 ice abundances vary little in the protostellar stage. Already, however, the 
low CH4 abundances do suggest significant ice destruction during the disk stage, through chemistry 
or selective desorption. Similarly, comets are known to be significantly depleted in nitrogen, and 
the small number of comets with directly measured NH3/H2O ratios reveal a loss of NH3 compared 
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to low mass protostars (JKawakita &: Mummall201ll . and references therein) - suggesting that the 



alteration of cometary ices is not specifically related to carbon. For a more complete cometary ice 
sample, see upcoming review by Mumma & Charnley (2011, ARAA review, in press). 



5. Conclusions 

Large samples of ice sources spanning diverse environments, evolutionary stages and luminosities 
are necessary both to determine 'typical' ice abundances and to identify how ice processes depend 
on their environment. The main findings from combining the c2d and other ice surveys carried out 
in a homogenous way are listed below: 

1. The CO and CO2 abundance medians relative to H2O are both ~29% toward low-mass YSOs, 
while CH4, NH3 and CH3OH are more than a factor of five less abundant. 

2. C02:H20, CH4 and NH3 abundances vary little with respect to H2O, suggesting co-formation. 
In contrast, C0:H20, OCN~, C02:C0 and CH3OH vary by factors 2-3 (lower to upper 
quartile) with respect to H2O, indicative of a separate formation pathway from H2O ice. Pure 
CO and CO2 ice abundances vary even more, consistent with their sensitivity to protostellar 
heating. 

3. Ice abundances toward background stars and protostars are similar except for a lack of features 
associated with ice heating, such as pure CO2 ice, toward background stars. In particular, 
there is no evidence for a different range of CH3OH and CO2 ice abundances. 

4. Protostellar ice abundances in nearby star-forming regions do not vary significantly between 
different clouds, except for a significantly lower pure CO abundances toward the protostars 
in Taurus. 

5. Compared to low-mass YSOs, ice abundances toward high-mass YSOs are characterized by 
low levels of CO and CO2 ice, which can be explained by ice heating. Low-mass protostars in 
the highly irradiated region IC 1396A have typical low-mass protostellar CO2 ice abundances, 
confirming that the original ice conditions during low- and high-mass star formation are 
similar. 

6. Correlation studies within the low-mass protostellar sample show a close association between 
CO, C02:C0, C0:H20 and the XCN band, supporting the latter's identification with OCN". 

7. Combining the above information, ice formation can generally be divided into three stages: an 
early phase driven by atomic hydrogenation reactions in clouds, which are fast compared to 
cloud core formation time scales of ~10^ years; a later CO-freeze-out dependent ice formation 
phase which takes place in the pre-stellar phase; and finally the protostellar phase where 
thermal and possibly UV processing shape the ice content. 



20 



8. Toward both low- and high-mass protostars and toward molecular clouds, an average of 8- 
16% of the total C, O and N are bound up in ices. Toward the most ice-rich source (IRS 51) 
60-80% of the non-refractory C and O are in ices. 

9. There is evidence for more complex ices toward both low- and high-mass protostars, e.g. 
HCOOH, CH3CHO and/or C2H5OH, but high-resolution spectra toward more sources are 
required to confirm their presence and quantify their abundances. 
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Fig. 1.— Ice spectra t oward the protostar s W33A (10^ L©), HH46 (12 Lq) and L1014 IRS (0.09 
Lq) (JGibb et al.l l200d : iBoogert et al.ll2004l . Paper I). The 3 /um portions of the spectra have been 
binned to increase the S/N. 
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Fig. 2. — Median abundances with respect to H2O in the low- and high-mass protostehar samples 
(^median) derived using survival analysis that includes upper limits (filled bars) and median detected 
abundances (outlined bars) toward low- and high- mass protostars. 
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Fig. 3. — Histograms of the CO (blue line-filled) and CO2 (black contours) abundances toward 
low-mass protostars with respect to the median abundances. The left panel shows the skewed 
abundance distributions of both molecules and the right panel the more symmetric log-normalized 
histograms. 
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Fig. 4. — Histograms of total ice abundances toward low-mass protostars (red) and high-mass 
protostars (black). For low-mass NH3 abundances, the solid filled histograms are from Paper IV 
using a silicate template to extract NH3 and the contoured histograms are from Paper IV and Table 
1 when using a polynomial to fit the local continuum around the NII3 feature. 
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Fig. 5. — Histograms of CO and CO2 components, otherwise as in Fig. 4. 
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Fig. 6. — Histograms of 5-7 fim components, otherwise as in Fig. 4. 
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Fig. 7. — Histograms of CO2 and the C1-C5 components (including upper limits) toward low-mass 
protostars (red contours), high-mass protostars (black contours) and the four IC 1396A sources 
(red solid). 
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Fig. 8. — Histograms of CO2 and CH3OH toward low-mass protostars (red fill), our background 
star sample (black line- fill) and Taurus background stars (blue contours). 



33 



0.8 



Q g _ CrM 2422 



o 

K 0.4 



0.2 - 



0.0 







A COgiCO X 5 
O CO 

^ 2175 cm"^ X 40 

+ C02:H20 




IRS44 



RS 46 



5.0x10^ 1.0x10^ 

Distance / AU 



1.5x10^ 
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Fig. 10. — Statistically significant correlations between ice abundances w. r. t. H2O and the 
H2O ice column density for low-mass protostars, with upper limits plotted with thin symbols. The 
significance of the correlations was measured using Spearman's rank correlation test, which makes 
no assumptions about the type of correlation, while the dotted line shows the best log-log fit to the 
data to guide the eye. 
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Fig. 11. — Panels a-d show correlations between XCN and CO ice components for low- mass pro- 
tostars (red diamonds) with high-mass protstellar abundances shown for comparison (black stars). 
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Table 1. Ice column densities, abundances and optical depths of the 5-7 //ni complex 

components. 



Source 


N(H20) 
1017 cm-2 


[NHs] 
% 


[CH3OH] 

% 


TCI (5.84 ^m) 


TC2 (6.18 ^lm) 


Tcs, (6.76 /im) 


TC4 (6.94 /im) 


TC5 (broad) 


WL 12 


22.1 ± 3.0 


<3.8 


<4.5 


0.012±0.003 


0.000±0.002 


0.049±0.004 


0.136±0.003 


0.014±0.071 


WL 6 


41.7 ± 6.0 


2.9 ± 0.4 


<2.1 


0.002±0.007 


0.000±0.006 


0.137±0.007 


0.087±0.006 


0.030±0.045 


IRS 42 


19.5 ± 2.0 


<2.1 


11.9 ± 1.1 


- 


- 


- 


- 


- 


IRS 43 


31.5 ±4.0 


- 


- 


0.059±0.004 


0.082±0.003 


0.179±0.005 


0.105±0.004 


0.066±0.051 


IRS 44 


34.0 ± 4.0 


3.7 ± 0.4 


<1.6 


0.080±0.005 


0.089±0.004 


0.135±0.005 


0.201±0.004 


0.004±0.086 


Elias 32 


17.9 ± 2.6 


<5.2 


12.4 ± 1.9 


0.021±0.005 


0.000±0.004 


0.050±0.009 


0.075±0.007 


0.058±0.025 


IRS 46 


12.8 ± 2.0 


5.1 ± 0.9 


<4.1 


0.018±0.004 


0.002±0.004 


0.076±0.005 


0.066±0.004 


0.000±0.023 


VSSG17 


17.0 ± 2.5 


<3.1 


6.9 ± 2.4 


0.058±0.002 


0.056±0.002 


0.042±0.005 


0.061±0.004 


0.011±0.016 


IRS 51 


22.1 ± 3.0 


2.4 ± 0.3 


11.7 ± 0.9 


0.042±0.003 


0.036±0.002 


0.074±0.003 


0.060±0.002 


0.000±0.012 


IRS 63 


20.4 ± 3.0 


5.7 ± 1.3 


<1.8 


0.000±0.003 


0.023±0.002 


0.048±0.003 


0.060±0.003 


0.044±0.039 
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Table 2. Typical ice abundances toward low- and high-mass protostars and background stars. 



Ice feature 


Low-mass 


High-mass 


background 


H2O 


100 


100 


100 


CO 


29 


13 


31 


CO2 


29 


13 


38 


CH3OH 


3 


4 


4 


NH3 


5 


5 




CH4 


5 


2 




XCN^ 


0.3 


0.6 





^XCN implies a band component securely identified with OCN and the nearby 2175 cm ^ component, which this study 

suggests is also due to OCN^. 
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Table 3. Abundance medians^ and lower and upper quartile values of ices and individual ice 

components with respect to H2O ice 



Ice feature 


Low-mass 


High-mass 


background 


H2O 




100 


100 


100 


CO 




3861 (29) 


13f 


31 


CO2 




29ii 


13i 


38M 


CH3OH 




7P (3) 


&f (4) 


810 (4) 


NH3 




04 


I6I0' (5) 






CH4 




5l 


4! (2) 






XCN 




0.6g:i (0.3) 


O.SJ'I (0.6) 






pure CO 




2lf 


3i 






C0:H20 




1319 


lop 






C0:C02 




2? 


1.3(V! (0.3) 






pure CO2 




2^3 


2? 






C02:H20 




20?i 


9f 


24 


C02:C0 




5l 


51 


6 


CO2 shoulder 


o.sji 


1? 




OCN- 




0.40:1 (0.2) 


o.eji 




2175 cm-i 




O.Sgi (0.2) 


101 




CI (HCOOH + H2CO)'' 


2.5f:i 


i2.8 
^■-^2.0 


2.8i:t 


C2 (HCOO- 


+NH3)'' 


1-9?:! (1-1) 


1 31-6 


2 ^3-3 


C3 (NH+ + 


CH30H)b 


4-3 i-l 


4.3i:l 


o 74.7 

•^■'3.4 


C4 (NH+)b 




2.2.1-1 


4-3i:§ 


2 l2-8 


C5 (warm H2O + anions)'' 


i.5g:i (0-9) 


3.36:1 (1.4) 





^'Values in parentheses include upper 
Since no single carrier the reported 



limits in the median calculation using survival analysis, 
number is peak optical depth / (A^HjO x lO^'^" X 100). 
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Table 4. Differences in CO and CO2 shoulder medians toward different star forming associations 





Pure CO 


CO2 


1 shoulder 


Ophiuchus 


32 




0.9 


Serpens 


38 




1.8 


CrA 


29 




1.8 


Perseus 


21 




0.8 


Taurus 


7 




0.5 


Other 


5 




0.5 


Complete sample 


21 




0.8 
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A. Ophiuchus NH3 and CH3OH spectra 

Figure HM shows the acquired Spitzer spectra for the nine young stellar objects in Ophiuchus 
that were added to the c2d sample. The spectra are dominated by a broad silicate feature which is 
removed by fitting a polynomial to regions free of molecular emission and absorption. We used the 
same fitting regions as in Paper IV to determine the ice band depths and shapes and then shifted 
the regions by ~0.1/im to determine the sensitivity of the results to the chosen local continuum. 

The resulting optical depth spectra are shown in Fig. [151 CH3OH ice is detected toward Elias 
32, IRS 42, IRS 51 and VSSG 17 and NH3 ice toward IRS 44, IRS 46, IRS 51, IRS 63 and WL 6. 
The bands are fitted with Gaussians and the derived peak positions and full-width-half-maxima are 
listed in Tabled The uncertainties are dominated by the choice of local continuum. The NH3 peak 
positions are consistent with laboratory ice spectra, while the FWHM are generally too low. From 
Paper IV, we know that the NH3 FWHM may be underestimated when using a local continuum 
rather than a silicate template to extract optical depth spectra and this explains the discrepancy. 
The CH3OH peak positions and FWHM are both consistent with laboratory measurements, with 
the peak positions suggesting that CH3OH is present in a CO-rich ice (Paper IV). 



Table 5. The amount of C, O and N bound up in protostellar ices. 



Cicc / Ctotal Oicc / Ototal Nice / Ntotal Cice / Cvol Oicc / Ovol Nice / Nvol 



Low-mass median 


15% 


16% 


10% 


27% 


34% 


10% 


Low- mass max 


46% 


29% 


35% 


83% 


61% 


35% 


High-mass median 


8% 


12% 


12% 


14% 


25% 


12% 


High-mass max 


18% 


17% 


22% 


32% 


36% 


22% 
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Wavelength {/J.ra) 



Fig. 14. — Spitzer spectra between 8.0 and 10.5 yum showing the sihcate feature and the superim- 
posed bands at 9 and 9.7 /im ascribed to NH3 and CH3OH ice, respectively. The red regions were 
used to fit a local continuum (blue dotted line). 
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Fig. 15. — Optical depth spectra of NH3 and CH3OH ice features toward the Ophiuchus sources 
obtained after subtracting the sihcate feature. The red hne shows the Gaussian fits to the detected 
bands. 
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Table 6: FWHM and peak positions of detected NH3 and CH3OH bands compared to laboratory 
values. 



Source 


Peak position NII3 


FWHM NH3 


Peak position CH3OH 


FWHM CH3OH 




fim 


cm ^ 


cm~^ 


^m 


cm ^ 


cm~^ 


Elias 32 








9.68±0.01 


1033±1 


36 ± 3 


IRS 42 








9.64±0.01 


1037±1 


47 ± 3 


IRS 44 


9.04±0.01 


1107±1 


28 ± 2 








IRS 46 


9.00±0.01 


lllOil 


27 ± 3 








IRS 51 


9.02±0.01 


1109±1 


25 ± 1 


9.70±0.01 


1031±1 


38 ± 2 


IRS 63 


9.07±0.02 


1102±2 


32 ± 5 








VSSG 17 








9.71±0.02 


1030±3 


45 ± 10 


WL 12 














WL 6 


9.04±0.01 


1106±1 


31 ± 2 








Laboratory'' 


8.85-9.41 


1062-1130 


53-137 


9.67-9.83 


1017-1034 


22-39 




) 
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B. Carriers of the 7.25 fiui ice feature 

Paper I showed that it is difficult to fit the 7.25 /im with spectra of HCOOH, its most com- 
monly assigned carrier. Laboratory spectra of some tertiary mixtures come closest, but suffer from 
uncertainties in baseline determinations (see below). 

Figure 1161 shows the spectr a of pure HCOOH ice and an H^O:HCOO H 4:1 ice mixture (a t 
15 K from bisschop et aLllJonTl ^). pure CH3CH2OH and CH3CHO spectra (|6berg et al.1 bonQbl ) . 
and observed spectra toward W33A, NGC7538 IRS9 and Bl-b. W33A and Bl-b clearly have 
absorption features at 7.25 and 7.4 //m, while NGC7538 IRS9 does not. The comparison shows 
that CH3CH2OH is a plausible carrier for the 7.25 ^m feature, commonly assigned to HCOOH. The 
lack of the 7.25 //m to ward NGC7538 IR S9 is consistent with the upper limit on CH3CH2OH ice 
from the 3 //m feature (JBoudin et al.lll998l ). while no such limits are published for the other sources. 
Figure [16] also shows that the observed 7.40 /xm feature is likely due to CH3CHO. Thus there is 
also little evidence for HCOO~ - both have been proposed previously as carriers (e.g. ISchutte et al 



I999I : iGibb et al.ll2004l . Paper I) 



These tentative identifications may in the future provide an opportunity to observe a complex 
ice chemistry in situ, but the Spitzer-IRS spectra are too low in spectral resolution to securely 
distinguish between different carriers, especially since the HCOOH spectral feature depends on the 
ice matrix and has been reported to become as narrow as the CH3CH2OH band i n mixtures with 



H2O and CH3OH, but this result depends crucially on the choice of local baseline (JBisschop et al 



20071 ). Until the HCOOH baseline issue has been resolved and higher resolution spectra exist for 
more sources the HCOOH and HCOO~ abundances should not be determined from the 7.25 and 
7.42 //m ice features. This does not imply that HCOOH and HCOO" are not present in interstellar 
ices, only that we do not currently have the tools to quantify their abundances. 
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Fig. 16.— ISO spectra at 7-8 ^um for W33A, NGC7538 IRS9 and Bl-b, following subtraction 
of a local spline continuum, plotted together with laboratory spectra of pure HCOOH, CH3CHO 
and CH3CH2OH ices. The dashed lines mark the 7.25 and 7.40 fim features usually assigned to 
HCOOH and HCOO~. The feature at 7.67 /um is due to CH4 ice. The baselines were obtained by 
fitting a local spline continuum to 6.9, 7.16, 7.33, 7.77 and 7.85 /im. 
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C. Histograms of ice detections and detections+upper limits 




Fig. 17. — Histograms of total ice abundances toward low-mass protostars (red) and iiigh-mass 
protostars (black). Solid iiistograms mark detections and contours detections+upper limits. For 
low-mass NH3 the solids are from Paper IV using a silicate template to extract NH3 and the contour 
are from Paper IV and Table 1 when using a polynomial to fit the local continuum. 
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Fig. 18. — Histograms of CO and CO2 components, otherwise as in Fig. [T71 
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Fig. 19. — Histograms of 5-7 fxni components, otherwise as in Fig. [T71 



